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ABSTRACT 

We analyze spatially resolved spectroscopic observations of the Eta Carinae binary system obtained 
with HST/STIS. Eta Car is enshrouded by the dusty Homunculus nebula, which scatters light emitted 
by the central binary and provides a unique opportunity to study a massive binary system from 
different vantage points. We investigate the latitudinal and azimuthal dependence of Ha line profiles 
caused by the presence of a wind- wind collision (WWC) cavity created by the companion star. Using 
two-dimensional radiative transfer models, we find that the wind cavity can qualitatively explain the 
observed line profiles around apastron. Regions of the Homunculus which scatter light that propagated 
through the WWC cavity show weaker or no Ha absorption. Regions scattering light that propagated 
through a significant portion of the primary wind show stronger P Cygni absorption. Our models 
overestimate the Ha absorption formed in the primary wind, which we attribute to photoionization 
by the companion, not presently included in the models. We can qualitatively explain the latitudinal 
changes that occur during periastron, shedding light on the nature of Eta Car's spectroscopic events. 
Our models support the idea that during the brief period of time around periastron when the primary 
wind flows unimpeded toward the observer, Ha absorption occurs in directions toward the central 
object and Homunculus SE pole, but not toward equatorial regions close to the Weigelt blobs. We 
suggest that observed latitudinal and azimuthal variations are dominated by the companion star via 
the WWC cavity, rather than by rapid rotation of the primary star. 

Keywords: stars: atmospheres — stars: mass- loss — stars: variables: general — supergiants — stars: 
individual (Eta Carinae) — binaries: general 



1. INTRODUCTION 

Eta Carinae is a rarity among the rarities that consti- 
tute massive stars in the short, unstable Luminous Blue 
Variable (LBV) phase. The nature and fate of the central 
object are still the subject of intense debates, although 
consensus seems to exist around an ecce ntric massive bi- 
nary scenario, as initially proposed by iDamineli et al.l 
(1997). This scenario is s upported by multi-wa velength 
observations from X-ravs (I Corcoran et al.l l2010) to radio 
wavelengths (|Duncan L White! I2003f ). 

The ultraviolet and optical spectra are strongly dom- 
inated by the dense wind of the primary star (hereafter 
77a) 5 hinting that most of the luminosity of the sys- 
tem (L* ~ 5x 1 Q 6 L(7), [ Davidson & Humphrevs 1997T) 
come s from t]a (Hil lier et al.l I2QQ1L 120061 : iGroh et all 
1201 2l hereafter G12). Recent spectroscopic analysis sug- 
gests that t]a has a mass- loss rate of 8.5 x lO _4 M yr _1 
and wind terminal velocity of ~ 420 krns" 1 (G12; see 
also lHillier et al.ll2Q0lLl2006h . The companion star (here- 
after t]b) has yet to be observed directly, with only in- 
direct constraints available on its temperature and lu- 
minosity (r ef f ~ 36,000 - 41, 000 K and 10 5 L© < 
L* < 10 6 L©; iMehner et al.l l2010h . and wind properties 
(voo ~ 3000 krns -1 and M - 1.4 x 10" 5 M^r -1 ; e.g., 
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iParkin et al1l20TTh . 

Eta Car has interested astronomers ever since its dra- 
matic brightness increase in the 1830s. This eruptive 
even t caused the ejectio n of tens of solar masses of mate- 
rial (jSmith et al.ll20Q3bh . enshrouding the central object 
within a dusty nebula (the Homunculus). The expand- 
ing dust residing in the Homunculus scatters light emit- 
ted by the central so urce. The 3D determ i nation of the 
Homunculus' shape (jDavidson et al.l 120011 : ISmithl l2Q06h 
allows one to relate the scattered light at a given posi- 
tion in the Homunculus to a line-of-sight t hat views the 
stellar system from a certain latitude (jSmith et a l. 2003a, 
hereafter S03) and azimuth. In the context of a binary, 
this implies that different positions in the Homunculus 
see the binary as if viewed from different orbital inclina- 
tions i and longitudes of periastron lj (Fig. [I]). 

Spatially resolved observations of Eta Car show that, 
during most phases of the orbital cycle, H and Fe II 
spectral lines are latitude dependent. The P Cygni ab- 
sorption is stronger in spectra taken around the pole 
of the Homun culus SE lobe than in the spectra of the 
central object (iHpexjTA llen 1992; Smith et al. 2003a : 
iWeis et all [20051 : IStahl etall 120051 : IMehner et all [2012 . 
K. Nielsen 2009, priv. comm.). The P Cygni absorption 
seems to extend to higher velocities in spectra reflected 
off the Homunculus SE pole, which, combined with the 
stronger absorption, has been interpreted as evidence for 
a denser, fa ster polar wind generated by the rapid rota- 
tion of t]a (|S03l ). The latitudinal dependence of H and 
Fe II seems to be less pronounced around periastron, sug- 
gesting that 77A 's wind becomes more spherical (jS03l ). 

However, H and Fe II lines, which are formed in the 
wind of 77a, are extremely affected by the presence of 
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Figure 1. Offset positions discussed in this Letter. The angles 
i and uo correspond approximately to the equivalent orbital incli- 
nation (i. e. latitude) and longitude of periastron (i. e. azimuth) 
that light scattered off different regions of the Homunculus see the 
central binary system. The background image shows a ~ 20" box 
around the central source in Eta Car (credit: Nathan Smith (U. 
Arizona), NASA). 

a cavity created by t?b (G12). This cavity arises nat- 
urally in hydro dynamical models as a consequence of 
the ca rving of the wind of 7?a by the fast, thin wind 
of 7?b (|Pittard fc Corcoranll2QQ2[ ). Taking into account 
the presence of the low-density cavity, two-dimensional 
(2D) radiative transfer models of Eta Car show that H 
and F e II lines become latitude- and azimuth-dependent 
([G12h . This occurs because of the varying amount of pri- 
mary wind material towards different latitudes and az- 
imuths. These models reproduce the absence of P Cygni 
absorption in Ha and Fe II lines in spectra seen in the di- 
rect view to the central object (which is viewed through 
the rarified cavity), without evoking rapid rotation of t]a- 
Can we still consider t]a as a rapid rotator and proto- 
type of massive stars with dense polar winds? The goal 
of this Letter is to investigate the origin of the observed 
latitudinal and azimuthal dependencies of t?a's wind and, 
in particular, whether they can be explained by the pres- 
ence of a WWC cavity in the wind of t]a- 

2. OBSERVATIONS AND MODELING 

To illustrate our findings, we focus on archival 
HST/STIS observations obtained around apastron (2000 
Mar 20, orbital phase <j> = 10.40fl 52" x 0'/2 aper- 
ture) and periastron (1998 Mar 20, <j> = 10.046, 52" x 
O'/l aperture). Both datasets were acquired with the 
G750M grating, covering the region around Ha with 
R ~ 8000. The spectra were extracted using custom 
IDL routines using an aperture of 6 half-pixels (0 / /152) on 
the central object and close to the Weigelt blobs (WBs; 
IWeieelt fc Ebersbergerifl98l located 0'/15-0'/30 NW of 
the central source), and 12 half-pixels elsewhere. The 
spectra were corrected for the velocity shift introduced 
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by the dust scatt ering using narrow f orbidden lines as 
proxies, following iMehner et al.l (|2012f ). This technique 
could not be applied to the offsets at 5"NE and 4 // SW, 
and the profiles were aligned on the red s ide of Ha, where 
the influence from t]b is minimized (|G12[ ) . Note that this 
does not affect our conclusions. We refer the reader to 
e. g. [§03 for details about the observations and data re- 
duction. 

The aforementioned observations are analyzed using 
the 2D radiative transfer models of Eta Car from G12. 
The 2D models take into account t]a and the presence 
of a rarified WWC cavity and a dense wind-wind inter- 
acting region, corresponding to the post-shocked primary 
wind. We refer the reader to lG12l for an extensive discus- 
sion about the impact of the wind- wind collision (WWC) 
cavity on the spectrum, and to their table 2 for the full 
model parameters. We assume an orbit al orientation 
with i = 138°, uo = 270°, PA Z = 312° (jMadura et al.l 
2012), and half-opening angle of the cavity of a = 57° 
(|G12[ ) . This corresponds to the direct view of the central 
source. Since offset positions in the Homunculus see the 
binary system under different vantage points, 2D mod- 
els with different i and oj were computed for each offset 
analyzed here. 

For visualization purposes, we also present 3D render- 
ings of Smoothed Particle Hydrodynamics (SPH) simu- 
lations of t he Eta Car binary syste m, similar to those 
presented in lMadura fc Grohl (|2012f ). 

3. LATITUDINAL CHANGES IN LINE PROFILES DUE TO 
THE WWC CAVITY 

We first investigate the strong latitudinal changes in 
line profiles around apastron. In the WWC cavity sce- 
nario, the latitudinal variations of line profiles depend on 
the amo unt o f primary wind material towards a certain 
latitude (|G12f ). This is regulated mainly by the size of 
the line formation region compared to the distance of the 
WWC apex to t?a, meaning that certain lines (Ha, Fe II) 
are more affected by the cavity than others (e.g., higher 
Balmer lines). For brevity, we present here results for 
Ha, which is formed at a distance of ~ 10 — 100 AU from 
7?a- Similar conclusions would be obtained for Hf3 and 
Fe II lines. 

Figure [2fc,d presents the observed and 2D model spec- 
tra of Eta Car for different viewing angles, corresponding 
to a direct view of the stellar system (black line), a view 
from the Homunculus SE pole (light blue), and from the 
Homunculus equator towards the WBs (red). Our 2D 
models qualitatively reproduce the observed variation of 
Ha with latitude, suggesting that the latitudinal changes 
can be explained by the presence of the WWC cavity. 
For viewing angles corresponding to line-of-sights that 
view t]a through the rarified cavity, there is little or no P 
Cygni absorption because of the reduced density within 
the cavity. This drastically reduces the n = 2 level popu- 
lation of H, reducing the Ha absorption. This is the case 
for the direct view to the central source or to the Ho- 
munculus equator towards the WBs (Fig. [2^). There is 
also a fair quantitative agreement between observations 
and 2D models at these two positions. 

For line-of-sights that cross the undisturbed wind of 
t?a, such as for the spectra reflected on the Homunculus 
SE pole, the population of the n = 2 energy level of H is 
unchanged in our model compared to a spherically sym- 
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Figure 2. a: 3D rendering of hydrodynamical simula- 
tions showing the extended, optically-thick primary wind 
(cyan/green/yellow/red) and the WWC zone (dark blue) at or- 
bital phase <j> = 10.409. Material in the primary wind is color 
coded to radius, i.e., red corresponds to material at larger radii 
(~ 75 — 150 AU) from the central stars, while green indicates ma- 
terial near the apex of the WWC region (< 15 AU from 77 a for 
these phases). Material below the orbital plane has been removed, 
and the 3D rendering slightly tilted, to allow better visualization of 
the dynamics of the inner regions. The physical scale corresponds 
to ~ ±150 AU centered on the center of mass of the binary sys- 
tem, b: Similar to (a), but for the 2D radiative transfer models, c: 
Observed spectra obtained at <f> = 10.409 at the Homunculus equa- 
tor (red), direct view at 138° (black), and Homunculus SE pole 
(blue). To aid the comparison between model and observations, 
the vertical gray dashed line corresponds to v = —500 kms -1 . d: 
Model spectra corresponding to these viewing angles. Note that 
the narrow emission seen in c around —40 kms -1 is due to nebular 
emission from the Weigelt blobs, and not included in the modeling. 

metric case. This results in stronger absorption of the 
stellar continuum in spectra scattered off the Homuncu- 
lus SE pole than in spectra of the direct view to the 
central source, in agreement with the observations. How- 
ever, our 2D model overestimates the amount of P Cygni 
absorption seen in spectra scattered off the Homunculus 
SE pole, and underestimates the maximum velocity of P 
Cygni absorption in Ha by ~ 70 kms -1 (blue lines in 
Fig. Efc,d). 

Our favored explanation for these discrepancies is that 
the polar r egions of 77 a 's wind could be mild ly photoion- 
ized by rj B (jMadura et al.ll2Q12t lKruipll2011h . This effect 
is not included in our 2D models and would reduce the 
population of the n = 2 energy level of H and thus the 
amount of Ha absorption. In addition, photoionization 
by t]b would likely change the driving of the wind of 
t?a, perhaps allowing it to reach velocities slightly higher 
than regions of t^a's wind not affected by t^b- 

Alternatively, a wrong choice of stellar and wind pa- 
rameters could also be the culprit, e. g. a lower M over 
the poles would be needed to produce less absorption 
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Figure 3. Similar to Fig. [2] but for azimuthal variations. 

in Ha, but the fit to the UV and optical emission lines 
would be significantly worse (|G12l ). In addition, we can- 
not rule out that the observations do not correspond to 
the pure polar spectrum, with the intrinsic P Cygni ab- 
sorption being diluted by continuum and Ha emission 
from the ejecta in the inner arcsecond. 

Lastly, we cannot preclude that the wind i s sli ghtly 
faster at the pole because of rapid rotation (S03), al- 
though our 2D cavity models suggest that Voo over the 
pole would be much smaller (~ 500 kms -1 ) than previ- 
ous observational estimates. 

4. AZIMUTHAL VARIATIONS DUE TO THE WWC CAVITY 

We turn our attention to azimuthal variations observed 
in scattered light off the Homunculus around apastron. 
These would not be expected in a single rapid rotator 
scenaricQ, while a WWC cavity intrinsically produces az- 
imuthal density variations (Fig. [3^). Therefore, this cru- 
cial difference could allow one to distinguish between the 
two scenarios. 

Figure [3fc presents the Ha line profiles scattered off the 
Homunculus equator, around the WBs (red), 5"NE (or- 
ange), and 4"SW (purple). These positions probe differ- 
ent azimuths, as if the system were viewed from i = 90°, 
but different uj. The corresponding synthetic Ha line pro- 
files from our 2D models are shown in Fig. [3]i. Our 2D 
model qualitatively explains the variations in strength of 
the P Cygni absorption as a function of azimuth for an 
orbital orientation with uj ~ 270°, i ~ 138°, and a = 57°. 
Lower values of uj would be allowed for larger a. 

The changes in the P Cygni absorption occur because 

7 A strong misalignment of the rotation and Homunculus polar 
axes could cause azimuthal variations in scattered light off the Ho- 
munculus. However, it would be challenging to produce P Cygni 
absorption both at the Homunculus pole and equator, as observed. 
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of the variation of primary wind material at different az- 
imuths in the equatorial regions. Line-of-sights that view 
the system down the WWC cavity have little primary 
wind material between the observer and t^a- This causes 
less absorption in H<x This is the case for the spec- 
tra scattered off the WBs (red) and at 4"SW (purple). 
Conversely, the spectrum scattered off the Homunculus 
equator at 5"NE probes essentially high densities corre- 
sponding to the unmodified primary wind. The quanti- 
tative agreement of our 2D models is satisfactory at WB 
and 4 // SW, but overestimates the amount of absorption 
seen at 5"NE. This discrepancy resembles that seen in 
the Homunculus polar spectrum and could also be caused 
by neglecting photoionization from t]b in our 2D models 
(Sect. E). 

5. VARIATIONS AROUND PERIASTRON AND THE 
NATURE OF THE SPECTROSCOPIC EVENTS 

Because Eta Car is a highly eccentric binary system, 
the carving of the primary wind by t]b changes as a func- 
tion of <j). Significant variations are expected around pe- 
riastron, when the cavity gets closer to t]a and becomes 
distorted due to the high orbital velo cities of t]b (e. g. 
IQkazaki et aHl2008t [Parkin et al.ll20TTh . 

Figure |4^ shows a 3D rendering of the primary wind 
and WWC cavity at <p — 10.046. Its morphology and ori- 
entation with respect to the observer are greatly modified 
compared to apastron (Fig. [2^). The observed line pro- 
files are still latitude-dependent, with the key difference 
being the appearance of P Cygni absorption in the spec- 
trum of the direct view to the central object (Fig. Hh). 
Notice that the spectrum scattered off the WBs (red) 
does not show Ha absorption near periastron (jGull et al.l 
[20(1). 

The hydrodynamics and morphology of the WWC cav- 
ity around periastron are intrinsically 3D and, as such, 
our 2D modeling fails to fully reproduce the observed 
changes at all latitudes and azimuths. However, as im- 
portant insights could be obtained with a 2D model, we 
present here a first attempt to model the periastron spec- 
tra at different latitudes. 

Figure |4]i displays the synthetic Ha line profiles from 
our 2D model, assuming the density structure depicted in 
Fig. HJd. After periastron, for a short amount of time the 
wind of t]a is able to escape to regions previously occu- 
pied by the WWC cavity, up to a distance of ~ 60 AU at 
</>orb=l 0.046. This material increases the column density 
of neutral H, in particular for viewing angles with i cor- 
responding to the direct view. Our 2D models indicate 
that this increase in primary wind material is enough to 
cause Ha P Cygni absorption (Fig. 0]i). 

This portion of primary wind material, which causes an 
increase in column density, could be easily confused with 
the shell ejectio n scenario advocated by earlier studies 
(jDav idson 199jJ S03). However, the key difference is that 
no change in M, nor instability in t]a, are needed to eject 
the required portion of primary wind material in line- 
of-sight to t]a- It is simply the result of the dynamical 
interplay between the primary and sec ondary wind in this 
eccentric system (Ma dura et al.ll2012h . For lines formed 
in t]a , we find that the portion of primary wind that flows 
unimpeded for a brief period of time after periastron has 
a significant impact on the line profiles. 

The fits to spectra scattered at the Homunculus SE 
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Figure 4. Similar to Fig. [2] but for just after periastron 
(<£ orb =10.046). 

pole suffer from a problem similar to that seen at apas- 
tron, when the P Cygni absorption is overestimated by 
our 2D model. However, our model can qualitatively re- 
produce the velocities of the P Cygni absorption seen in 
polar and direct view spectra. In addition, the 2D models 
predict little variability in the polar spectrum between 
apastron and periastron, i n agreement with t he obser- 
vations (jStahl et al.ll2QQ5l : iMehner et a71l2012[ ). The P 
Cygni absorption at the equator is also overestimated by 
the 2D model. We attribute this to the model assump- 
tion breaking down at this latitude, since the density 
structure of the primary wind is severely modified by the 
WWC cavity (Fig. g]). Namely, the 3D SPH simulations 
predict three structures in line-of-sight to t]a when seen 
from the equator: a fossil outer cavity, a geometrically 
thin (~ 20 AU) portion of primary wind, and a newly- 
formed cavity by tjb- The Ha absorption arising in our 
model comes from the thin portion of primary wind ma- 
terial, but the assumed source function and opacities are 
unrealistic. This is because the newly formed cavity al- 
lows photons from t]a and t]b to ionize this thin portion 
of primary wind in the equatorial regions. Because of 
our extremely simplified assumptions here, it is manda- 
tory to include 3D and ionization effects to improve the 
modeling at periastron. 

6. IMPLICATIONS: A COMPANION THAT DOES IT ALL? 
We showed that the 2D models of Eta Carinae from 



IG12I qualitatively explain the latitudinal and azimuthal 
variations of Ha line profiles observed around apastron. 
We also attempted a first modeling of the variations 
seen during periastron, showing promising results and ex- 
plaining the increase in Ha P Cygni absorption observed 
in the direct view to the central object after periastron. 
Latitudinal variations in P Cygni absorption lines have 
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been interpreted as being caused by the rapid rot atio n 
of 77A, which would produce a dense polar wind (S03). 
We suggest that an otherwise spherical t]a wind carved 
by t]b can explain these variations. The fact that the 
presence of t]b affects the Ha and Fe II line profiles (|G12l ) , 
which have been used as diagnostics of rapid rotation, 
hampers the determination of how rapid t]a is spinning. 
Ultimately, this also raises the point of whether rapid 
rotation is really present. 

More importantly, our 2D models reproduce not only 
the latitudinal, but also the azimuthal variations ob- 
served in spectra scattered off the Homunculus equator. 
Azimuthal changes are in principle not expected from 
rapid rotation, and additional mechanisms would have 
to be evoked to explain the observations. Our cavity 
scenario naturally produces azimuthal variations, caused 
by different amounts of primary wind material at differ- 
ent azimuths, which is controlled by the opening angle 
of the WWC cavity. The results from this Letter cast se- 
rious doubts on the idea that t]a has been a rapid rotator 
in recent decades. 

Our 2D models can qualitatively explain the latitudinal 
dependenc e variations se e n in t he ground-based Ha line 
profiles of lMehner etHI (|2012l ) (their Fig. 8). Our 2D 
models also suggest that the spectrum reflected off the 
Homunculus pole should remain essentially unchanged 
during periastron, in excellent agreement with the ob- 
servations. The lack of significan t long-term variation 
in the polar (Me hner et a l. 2012) and equ atorial spec- 
trum reflected off the WBs (|Gull et al.ll2009D . both in line 
strengths and continuum flux, also would argue against 
dramatic long-term changes in stellar parameters of tja- 

Finally, one might wonder how Eta Car fits in the con- 
text of stellar evolution thro ugh the LBV phase . It might 
be that AG Car or HR Car (|Groh et al.ll2QQ9al lbh. which 
so far do not have detected companions, might be more 
appropriate prototypes for Galactic LBVs rather than 
Eta Car. Also, how does Eta Car fare in the context 
of the most massive stars, such as those observed in the 
core of the massive st ellar clusters NGC 3603 and R136 
(jCrowther et al.l l2010)? Is Eta Car an evolved version of 
those stars, at a stage when a Giant Eruption recently 
occurred? Or is binary evolution required to reach the 
stage where Eta Car is today? Specially if the presence of 
t]b turns out to be a necessary factor for the occurrence 
of the Giant Eruption, perhaps the main observational 
properties that make Eta Car so unique in the Galaxy are 
directly or indirectly linked to the presence of a binary 
companion. 

Our results provide the basic foundation to understand 
line profile variations scattered off the Homunculus, and 
to disentangle the long-term evolution of the system from 
phase- locked changes. The presence of a WWC cavity 
seems to be the dominant effect to understand Ha on the 
central source and WBs. We encourage future theoretical 



efforts to take into account the ionization from which 
is expected to play a role at selected orbital phases and 
latitudes/ azimuths, and to explain the behavior of other 
lines such as US and He I. 
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